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Abstract 

We summarize the top-quark mass measurements from the CDF and D0 experiments at 
Fermilab. We combine published Run-I (1992-1996) measurements with the most recent 
preliminary Run-II (2001-present) measurements using up to 1 fb" 1 of data. Taking 
correlated uncertainties properly into account the resulting preliminary world average 
mass of the top quark is M t = 171.4 ± 1.2(stat) ± 1.8(syst) GeV/c 2 , which corresponds to 
a total uncertainty of 2.1 GeV/c 2 . The top-quark mass is now known with a precision of 
1.2%. 
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1 Introduction 



The experiments CDF and D0, taking data at the Tevatron proton-antiproton collider located 
at the Fermi National Accelerator Laboratory, have made several direct experimental measure- 
ments of the top-quark pole mass, M t . The pioneering measurements were based on about 
100 pb" 1 of Run-I (1992-1996) data HH^EHHElinilZlIHlISl HOI HH H3 and include results from 
the it — > qq'bqq'b (all-j), the it — > ivqq'bb (1+j), and the it — > £ + vbl~Vb (di-1) decay channels 2 . 
Results using approximately 350 pb _1 of Run-II (2001-present) data have been recently pub- 
lished in the 1+j and di-1 channels f3J El 113 • The Run-II measurements summarized here 
are the most recent preliminary results in the 1+j and di-1 channels using 370 - 1000 pb" 1 of 
data and improved analysis techniques [THl H3 [TBI CH] • This combination also includes a Run-II 
measurement in the all-j channel [2U| and a new analysis technique in the 1+j channel [2JJ (lxy ) . 



This note reports the world average top-quark mass obtained by combining five published 
Run-I measurements [2J El El El UHl El with the most recent preliminary Run-II measurements 
from CDF [12111711201121] and D0 [13111]. The combination takes into account the statistical 
and systematic uncertainties and their correlations using the method of references [221 122] and 
supersedes previous combinations [23 [23 I2E] • The most precise individual measurements of 
Mt are now the preliminary measurements in the 1+j channel from Run II. These are 170.9 ± 
2.5 GeV/c 2 (CDF, |1&) and 170.3 + 4.5 GeV/c 2 (D0, [JS]). These have weights in the new M t 
combination of 62% and 19%, respectively. 

The input measurements and error categories used in the combination are detailed in Sec- 
tion [2] and El respectively. The correlations used in the combination are discussed in Section 0] 
and the resulting world average top-quark mass is given in Section El A summary and outlook 
are presented in Sectional 



2 Input Measurements 

For this combination eleven measurements of M t are used, five published Run-I results and 
six preliminary Run-II results. In general, the Run-I measurements all have relatively large 
statistical uncertainties and their systematic uncertainty is dominated by the total jet energy 
scale (JES) uncertainty. In Run-II both CDF and D0 take advantage of the larger it samples 
available and employ new analysis techniques to reduce both these uncertainties. In particular 
the JES is constrained using an in-situ calibration based on the invariant mass of W — > qq' 
decays in the 1+j channel. The Run-II D0 analysis in the 1+j channel constrains the response of 

2 Here i = e or /i. Decay channels with explicit tau lepton identification are presently under study and are 
not yet used for measurements of the top-quark mass. 
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light-quark jets using the in-situ W — > gg' decays. The JES determined in this manner is then 
propagated to their measurement in the di-1 channel. Residual JES uncertainties associated 
with 77— and pr-dependencies as well as uncertainties specific to the response of 6-jets are 
treated separately. Similarly, the Run-II CDF analysis in the 1+j channel also constrains the 
JES using the in-situ W — ► qq' decays. Small residual JES uncertainties arising from rj— and p?- 
dependencies and the modeling of 6-jets are included in separate error categories. The Run-II 
CDF di-1 and all-j measurements use a JES determined from external calibration samples. 
Some parts of the associated uncertainty are correlated with the Run-I JES uncertainty as 
noted below. 

In previous combinations the Run-II CDF 1+j analysis used the JES determined from the 
external calibration as an additional gaussian constraint. This required us to treat that mea- 
surement as two separate inputs in the combination in order to accurately account for all the 
JES correlations. This gaussian constraint is not used in the present analysis as it does not 
significantly improve the sensitivity. Thus we can treat this measurement as a single input in 
the same manner as all the other measurements. 

A new analysis technique from CDF is also included (lxy). This measurement uses the 
mean decay-length from B-tagged jets to determine the top-quark mass. While the statisti- 
cal sensitivity is not nearly as good as the more traditional methods, this technique has the 
advantage that since it uses only tracking information, it is almost entirely independent of 
JES uncertainties. Additionally, since it does not require a full event reconstruction, it can 
use a more inclusive sample of tt candidates (e.g. events with > 3 jets). As the statitistics 
of this sample continue to grow, this method could offer a nice cross-check of the top-quark 
mass that's largely independent of the dominant JES systematic uncertainty which plagues the 
other measurements. The statistical correlation between this measurement and the Run-II CDF 
1+j measurement is determined using Monte Carlo signal-plus-background psuedo-experiments 
which correctly account for the sample overlap and is found to be consistent with zero (to 
within < 1%) independent of the assumed top-quark mass. 

The inputs used in the combination are summarized in Table ^ with their uncertainties 
sub-divided into the categories described in the next Section. The correlations between the 
inputs are described in Section EJ 



3 Error Categories 

We employ the same error categories as used for the previous world average |26j . They have 
evolved to include a detailed breakdown of the various sources of uncertainty and aim to lump 
together sources of systematic uncertainty that share the same or similar origin. For example, 
the "Signal" category discussed below includes the uncertainties from ISR, FSR, and PDF - 
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Run-I published 


Run-II preliminary 






CDF 




D0 






CDF 




D0 






all-j 


1+j 


di-1 


1+j 


di-1 


1+j 


di-1 


all-j 


lxy 


1+j 


di-1 


Lumi (pb 1 ) 


110 


105 


110 


125 


125 


955 


1030 


1020 


695 


370 


370 


Result 


186.0 


176.1 


167.4 


180.1 168.4 


170.9 


164.5 


174.0 


183.9 


170.3 178.1 


iJES 


0.0 


0.0 


0.0 


0.0 


0.0 


1.4 


0.0 


0.0 


0.0 


0.0 


0.0 


aJES 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


1.5 


bJES 


0.6 


0.6 


0.8 


0.7 


0.7 


0.6 


0.6 


0.5 


0.0 


0.6 


1.4 


cJES 


3.0 


2.7 


2.6 


2.0 


2.0 


0.0 


2.8 


3.1 


0.0 


0.0 


0.0 


dJES 


0.3 


0.7 


0.6 


0.0 


0.0 


0.2 


1.6 


0.8 


0.0 


3.5 


3.8 


rJES 


4.0 


3.4 


2.7 


2.5 


1.1 


0.0 


1.3 


3.0 


0.3 


0.0 


0.0 


Signal 


1.8 


2.6 


2.8 


1.1 


1.8 


1.1 


0.9 


0.9 


1.4 


0.5 


1.7 


BG 


1.7 


1.3 


0.3 


1.0 


1.1 


0.2 


0.7 


0.7 


2.3 


0.5 


1.0 


Fit 


0.6 


0.0 


0.7 


0.6 


1.1 


0.4 


0.9 


0.0 


4.8 


0.5 


0.9 


MC 


0.8 


0.1 


0.6 


0.0 


0.0 


0.2 


0.9 


1.0 


0.7 


0.0 


0.0 


UN/MI 


0.0 


0.0 


0.0 


1.3 


1.3 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Syst. 


5.7 


5.3 


4.9 


3.9 


3.6 


1.9 


3.9 


4.7 


5.6 


3.8 


4.8 


Stat. 


10.0 


5.1 


10.3 


3.6 


12.3 


1.6 


3.9 


2.2 


14.8 


2.5 


6.7 


Total 


11.5 


7.3 


11.4 


5.3 


12.8 


2.5 


5.6 


5.2 


15.8 


4.5 


8.3 



Table 1: Summary of the measurements used to determine the world average M t . All numbers 
are in GeV/c 2 . The error categories and their correlations are described in the text. The 
total systematic uncertainty and the total uncertainty are obtained by adding the relevant 
contributions in quadrature. 

all of which affect the modeling of the tt signal. Additional categories have been added in 
order to accommodate specific types of correlations. For example, the jet energy scale (JES) 
uncertainty is sub-divided into several components in order to more accurately accommodate 
our best estimate of the relevant correlations. Each error category is discussed below. 



Statistical: The statistical uncertainty associated with the M t determination. 

iJES: That part of the JES uncertainty which originates from in-situ calibration procedures 
and is uncorrelated among the measurements. In the combination reported here it cor- 
responds to the statistical uncertainty associated with the JES determination using the 
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W — > qq' invariant mass in the CDF Run-II 1+j measurement. Residual JES uncertain- 
ties, which arise from effects not considered in the in-situ calibration, are included in 
other categories. 

a JES: That part of the JES uncertainty which originates from differences in detector e/h 
response between fe-jets and light- quark jets. It is specific to the D0 Run-II measurements 
and is taken to be uncorrelated with the D0 Run-I and CDF measurements. 

bJES: That part of the JES uncertainty which originates from uncertainties specific to the 
modeling of 6-jets and which is correlated across all measurements. For both CDF and 
D0 this includes uncertainties arising from variations in the semi-leptonic branching 
fraction, ^-fragmentation modeling, and differences in the color flow between 6-jets and 
light-quark jets. These were determined from Run-II studies but back-propagated to the 
Run-I measurements, whose rJES uncertainties (see below) were then corrected in order 
to keep the total JES uncertainty constant. 

cJES: That part of the JES uncertainty which originates from modeling uncertainties corre- 
lated across all measurements. Specifically it includes the modeling uncertainties associ- 
ated with light-quark fragmentation and out-of-cone corrections. 

dJES: That part of the JES uncertainty which originates from limitations in the calibra- 
tion data samples used and which is correlated between measurements within the same 
data-taking period (ie. Run I or Run II) but not between experiments. For CDF this 
corresponds to uncertainties associated with the 7/-dependent JES corrections which are 
estimated using di-jet data events. For D0 Run-II this corresponds to uncertainties as- 
sociated with the light-quark response as determined using the W — > qq' invariant mass 
in the 1+j channel and propagated to the di-1 channel. The residual 77-dependent and 
PT-dependent uncertainties for the D0 Run-II measurements are also included here since 
they are constrained using Run-II 7+jet data samples. 

rJES: The remaining part of the JES uncertainty which is correlated between all measurements 
of the same experiment independent of data-taking period, but is uncorrelated between 
experiments. This is dominated by uncertainties in the calorimeter response to light- 
quark jets. For CDF this also includes small uncertainties associated with the multiple 
interaction and underlying event corrections. 

Signal: The systematic uncertainty arising from uncertainties in the modeling of the tt signal 
which is correlated across all measurements. This includes uncertainties from variations 
in the ISR, FSR, and PDF descriptions used to generate the it Monte Carlo samples 
that calibrate each method. It also includes small uncertainties associated with biases 
associated with the identification of 6-jets. 

Background: The systematic uncertainty arising from uncertainties in modeling the domi- 
nant background sources and correlated across all measurements in the same channel. 
These include uncertainties on the background composition and shape. In particular 
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uncertainties associated with the modeling of the QCD multi-jet background (all-j and 
1+j), uncertainties associated with the modeling of the Drell-Yan background (di-1), and 
uncertainties associated with variations of the fragmentation scale used to model W+jets 
background (all channels) are included. 

Fit: The systematic uncertainty arising from any source specific to a particular fit method, 
including the finite Monte Carlo statistics available to calibrate each method. 

Monte Carlo: The systematic uncertainty associated with variations of the physics model 
used to calibrate the fit methods and correlated across all measurements. For CDF it 
includes variations observed when substituting PYTHIA [27J EH EH! (Run I and Run II) 
or ISAJET 30 j (Run I) for HERWIG jSH E2] when modeling the it signal. Similar 
variations are included for the D0 Run-I measurements. The D0 Run-II measurements 
use ALPGEN [33] to model the it signal and the variations considered are included in 
the Signal category above. 

UN/MI: This is specific to D0 and includes the uncertainty arising from uranium noise in the 
D0 calorimeter and from the multiple interaction corrections to the JES. For D0 Run-I 
these uncertainties were sizable, while for Run-II owing to the shorter integration time 
and in-situ JES determination, these uncertainties are negligible. 



These categories represent the current preliminary understanding of the various sources of 
uncertainty and their correlations. We expect these to evolve as we continue to probe each 
method's sensitivity to the various systematic sources with ever improving precision. Varia- 
tions in the assignment of uncertainties to the error categories, in the back-propagation of the 
bJES uncertainties to Run-I measurements, in the approximations made to symmetrize the 
uncertainties used in the combination, and in the assumed magnitude of the correlations all 
negligibly effect (<C 0.1GeV/c 2 ) the combined M t and total uncertainty. 



4 Correlations 

The following correlations are used when making the combination: 

• The uncertainties in the Statistical, Fit, and iJES categories are taken to be uncorrelated 
among the measurements. 

• The uncertainties in the aJES and dJES categories are taken to be 100% correlated among 
all Run-I and all Run-II measurements on the same experiment, but uncorrelated between 
Run I and Run II and uncorrelated between the experiments. 

• The uncertainties in the rJES and UN/MI categories are taken to be 100% correlated 
among all measurements on the same experiment. 
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Run-I published 


Run-II preliminary 


CDF 
1+j di-1 all-j 


D0 
1+j di-1 


CDF 

1+j di-1 all-j lxy 


D0 
1+j di-1 


CDF-I 1+j 
CDF-I di-1 
CDF-I all-j 


1.00 

0.29 1.00 

0.32 0.19 1.00 








D0-I 1+j 
D0-I di-1 


0.26 0.15 0.14 
0.11 0.08 0.07 


1.00 

0.16 1.00 






CDF-II 1+j 
CDF-II di-1 
CDF-II all-j 
CDF-II lxy 


0.19 0.13 0.08 
0.36 0.23 0.26 
0.56 0.33 0.42 
0.07 0.03 0.02 


0.14 0.07 
0.24 0.12 
0.28 0.12 
0.05 0.01 


1.00 

0.13 1.00 

0.12 0.55 1.00 

0.05 0.03 0.03 1.00 




D0-II 1+j 
D0-II di-1 


0.07 0.04 0.02 
0.09 0.07 0.04 


0.06 0.02 
0.07 0.05 


0.09 0.03 0.03 0.03 
0.13 0.07 0.05 0.02 


1.00 

0.44 1.00 



Table 2: The resulting matrix of total correlation coefficients used to determined the world 
average top quark mass. 

• The uncertainties in the Background category are taken to be 100% correlated among all 
measurements in the same channel. 



• The uncertainties in the bJES, cJES, Signal, and Generator categories are taken to be 
100% correlated among all measurements. 



Using the inputs from Tableland the correlations specified here, the resulting matrix of total 
correlation co-efficients is given in Tabled 

The measurements are combined using a program implementing a numerical y 2 minimiza- 
tion as well as the analytic BLUE method [221 EH] ■ The two methods used are mathematically 
equivalent, and are also equivalent to the method used in an older combination [31], and give 
identical results for the combination. In addition, the BLUE method yields the decomposition 
of the error on the average in terms of the error categories specified for the input measure- 
ments [23|. 
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Run-I published 


Run-II preliminary 


CDF 

1+j di-1 all-j 


D0 
1+j di-1 


CDF 

1+j di-1 all-j lxy 


D0 
1+j di-1 


Pull 

Weight [%] 


+0.68 -0.35 +1.29 
-3.1 -0.6 -0.3 


+ 1.80 -0.23 
+8.0 +0.6 


-0.34 -1.33 +0.56 +0.80 
+61.7 +4.8 +10.3 +0.9 


-0.26 +0.85 
+ 18.9 -1.1 



Table 3: The pull and weight for each of the inputs used to determine the world average mass 
of the top quark. See Reference [22] 101 a discussion of negative weights. 

5 Results 

The combined value for the top-quark mass is: 

M t = 171.4 ± 2.1 GeV/c 2 , (1) 

with a x 2 °f 10.6 for 10 degrees of freedom, which corresponds to a probability of 39% indicating 
good agreement among all the input measurements. The total uncertainty can be sub-divided 
into the contributions from the various error categories as: Statistical (±1.2), total JES (±1.4), 
Signal (±0.9), Background (±0.3), Fit (±0.3), Monte Carlo (±0.3), and UN/MI (±0.1), for a 
total Systematic (±1.7), where all numbers are in units of GeV/c 2 . The pull and weight for each 
of the inputs are listed in Table El The input measurements and the resulting world average 
mass of the top quark are summarized in Figure H 

The weights of many of the Run-I measurements are negative. In general, this situation 
can occur if the correlation between two measurements is larger than the ratio of their total 
uncertainties. This is indeed the case here. In these instances the less precise measurement 
will usually acquire a negative weight. While a weight of zero means that a particular input 
is effectively ignored in the combination, a negative weight means that it affects the resulting 
central value and helps reduce the total uncertainty. See reference j22] for further discussion of 
negative weights. 

Although the x 2 from the combination of all measurements indicates that there is good 
agreement among them, and no input has an anomalously large pull, it is still interesting to 
also fit for the top-quark mass in the all-j, 1+j, and di-1 channels separately. We use the same 
methodology, inputs, error categories, and correlations as described above, but fit for the three 
physical observables, M t all_ -', M t 1+ -', and M t dl_1 . The results of this combination are shown in 
Tableland have x 2 °f 8-4 for 8 degrees of freedom, which corresponds to a probability of 40%. 
These results differ from a naive combination, where only the measurements in a given channel 
contribute to the M t determination in that channel, since the combination here fully accounts 
for all correlations, including those which cross-correlate the different channels. Using the 
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Parameter 


Value (GeV/c 2 ) 


Correlations 


M all-j 

M t 1+j 
M? 1 - 1 


173.4 ±4.3 
171.3 ±2.2 
167.0 ±4.3 


1.00 

0.29 1.00 

0.46 0.37 1.00 



Table 4: Summary of the combination of the nine measurements by CDF and D0 in terms of 
three physical quantities, the mass of the top quark in the all-jets, lepton±jets, and di-lepton 
channel. 

results of Table |U we calculate the chi-squared consistency between any two channels, including 
all correlations, as ^{dil — Ij) = 1.2, x 2 {lj ~ allj) = 0.24, and ^{allj — dil) = 2.1. These 
correspond to chi-squared probabilities of 27%, 65%, and 16%, respectively, and indicate that 
the determinations of M t from the three channels are consistent with one another. 



6 Summary 

A preliminary combination of measurements of the mass of the top quark from the Teva- 
tron experiments CDF and D0 is presented. The combination includes five published Run-I 
measurements and six preliminary Run-II measurements. Taking into account the statistical 
and systematic uncertainties and their correlations, the preliminary world-average result is: 
Af t = 171.4 ±2.1 GeV/c 2 . 

The mass of the top quark is now known with an accuracy of 1.2%, limited by the systematic 
uncertainties, which are dominated by the jet energy scale uncertainty. This systematic is 
expected to improve as larger data sets are collected since new analysis techniques constrain 
the jet energy scale using in-situ W — > qq' decays. It can be reasonably expected that with the 
full Run-II data set the top-quark mass could be known to better than 1%. To reach this level of 
precision further work is required to determine more accurately the various correlations present, 
and to understand more precisely the 6-jet modeling, Signal, and Background uncertainties 
which may limit the sensitivity at larger data sets. Limitations of the Monte Carlo generators 
used to calibrate each fit method may also become important as the precision reaches the 1% 
level and will warrant further study in the future. 



9 



Mass of the Top Quark (*Preliminary) 
Measurement 



CDF-I di-l 
D0-I di-l 
CDF-II di-l* 
D0-II di-l* 
CDF-I l+j 
D0-I l+j 
CDF-II l+j* 
D0-II l+j* 
CDF-I all-j 
CDF-II all-j v 
CDF-II Ixy* 



Tevatron Run- 



M top [GeV/c*] 



167.4 ± 11.4 

1 68.4 ± 12.8 
164.5± 5.6 
178.1 ± 8.3 
176.1 ± 7.3 
180.1 ± 5.3 
170.9± 2.5 
170.3± 4.5 

186.0 ± 11.5 
174.0± 5.2 

■183.9 ± 15.8 



X /dof = 10.6/10 



171.4 ± 2.1 



150 170 

M top [GeV/c 2 ] 



190 



Figure 1: A summary of the input measurements and resulting world average mass of the top 
quark. 
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